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Abstract. We discuss the basic features of the propagation of Ultra High Energy Cos- 
mic Rays in astrophysical backgrounds, comparing two alternative computation schemes 
to compute the expected fluxes. We also discuss the issue of the transition among galac- 
tic and extra-galactic cosmic rays using theoretical results on fluxes to compare different 
models. 

1 Introduction 

Ultra High Energy Cosmic Rays (UHECR) are the most energetic particles observed in nature with 
(detected) energies up to 3 H- 5 x 10 20 eV. The experimental observations of UHECR are performed 
nowadays by the Auger observatory in Argentina [ 1 ] and by the HiRes |2] and Telescope Array (TA) 
observatories in the USA. The observation of these particles poses many interesting questions 
mainly on their nature and origin. 

The study of the propagation of UHECR through astrophysical backgrounds can be very useful 
to interpret the observations and, may be, to have important insights on the possible sources of these 
particles. The propagation of UHECR from the sources to the observer is mainly conditioned by the in- 
teractions with the intervening astrophysical backgrounds such as the Cosmic Microwave Background 
(CMB) and the Extragalactic Background Light (EBL). While the propagation of protons is condi- 
tioned only by the CMB radiation field [4j, nuclei propagation is also affected by the EBL radiation 
0. 

Several propagation dependent features in the spectrum can be directly linked to the chemical com- 
position of UHECR and/or to the distribution of their sources |5|. Among such features particularly 
important is the Greisin, Zatsepin and Kuzmin (GZK) suppression of the flux, an abrupt depletion 
of the observed proton spectrum due to the interaction of the UHE protons with the CMB radiation 
field [6 |. The GZK suppression, as follows from the original papers, is referred to protons and it is 
due to the photo-pion production process on the CMB radiation field (p + jcub —> n + p). In the 
case of nuclei the expected flux also shows a suppression at the highest energies that, depending on 
the nuclei specie, is due to the photo-disintegration process on the CMB and EBL radiation fields 
(A + jcmb,ebl — * (A - nN) + nN) In any case, the interaction processes between UHE particles 
and astrophysical backgrounds will condition the end of the CR spectrum at the highest energies and 
the high energy behavior of the flux can be used as a diagnostic tool for the chemical composition of 
the observed particles. Another important feature in the spectrum that can be directly linked with the 
nature of the primary particles and their origin (galactic/extra-galactic) is the pair-production dip 0. 
This feature is present only in the spectrum of UHE extragalactic protons and, as the GZK, is a direct 
consequence of the proton interaction with the CMB radiation field, in particular the dip brings a direct 
imprint of the pair production process p + jcmb — * p + e + + e~ suffered by protons in their interaction 
with CMB radiation. 
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The experimental observations show contradictory results mainly on the chemical composition of 
UHECR. While Auger favors a light (proton) composition at low energies and an heavy (nuclei) com- 
position at the highest energies J8], HiRes and TA show a proton dominated spectrum at all energies 
(S). A clear understanding of the chemical composition of UHECR is of paramount importance in 
the study of such particles in particular in the determination of their possible sources and to tag the 
transition among galactic and extra-galactic CR J9]. 

Apart from the observations on chemical composition, also the observations on fluxes show a 
certain level of disagreement among HiRes and TA on one side and Auger on the other side. HiRes 
and TA experiments show a flux with a clear observation of the protons GZK suppression and the 
pair-production dip, coherently with their composition observations 1 10 1. The situation changes if the 
Auger results are taken into account. The latest release of the Auger data on flux shows a spectrum not 
compatible with the pair production dip with an high energy suppression not compatible with the GZK 
suppression as expected for protons. Signaling a possible deviation from a proton dominated spectrum 
to an heavier composition at the highest energies, confirming with the flux behavior the observations 
on chemical composition 1 10 1. 

This puzzling situation, with different experiments favoring different scenarios, shows the impor- 
tance of a systematic study of the propagation of UHECR in astrophysical backgrounds. In the present 
paper we will present two alternative ways to study the propagation of UHECR in astrophysical back- 
grounds: through a Monte Carlo (MC) approach and through an analytic computation scheme based on 
the kinetic theory of particles propagation. We will than compare theoretical results with observations, 
considering in particular the issue of the transition from galactic to extra-galactic CR. 

The paper is organized as follows: in section|2]we will briefly review theoretical models to study 
UHECR propagation, in section [3] we will use the results of the propagation models to determine the 
features of the transition among galactic and extra-galactic CR, in this section we will use the results 
of lfT31l for the galactic component, finally we will conclude in section|4] 



2 Ultra High Energy Cosmic Rays Propagation 

The propagation of charged particles (protons or nuclei) with energies above 10 17 eV through astro- 
physical backgrounds can be suitably studied taking into account the main channels of interaction that, 
as already anticipated in the introduction, are: 

- protons UHE protons interact only with the CMB radiation field giving rise to the two processes 
of pair production and photo-pion production. 

- nuclei UHE nuclei interact with the CMB and EBL radiation fields, suffering the process of pair 
production, in which only CMB is relevant, and photo-disintegration, that involves both back- 
grounds. While the first process conserves the nuclear species, the second produces a change in 
the nuclear species, extracting nucleons from the nucleus. 

In the energy range E ^ 10 18 -r 10 19 eV the propagation of UHE particles is extended over cos- 
mological distances with a typical path length of the order of Gpc. Therefore we should also take into 
account the adiabatic energy losses suffered by particles because of the cosmological expansion of the 
Universe. 

In this section we will briefly recall two alternative computation schemes used to describe the 
propagation of UHECR: the kinetic approach [4] and the SimProp MC approach IfTTII . 

The kinetic approach is based on the hypothesis of continuum energy losses (CEL), through which 
particles interactions are treated as a continuum process that continuously depletes the particles energy. 
In the propagation through astrophysical backgrounds the interaction of particles are naturally affected 
by fluctuations, with a non-zero probability for a particle to travel without loosing energy. The effect 
of the fluctuations in the interaction of particles is taken into account in the MC approach while it is 
neglected in the case of the kinetic approach (CEL approximation). In the case of proton propagation 
the CEL approximation has a negligible effect on the pair-production process, while in the case of 
photo-pion production it gives a deviation only at the highest energies (E > 10 20 eV) of the order of 
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Fig. 1. Flux of iron and secondary nuclei at z = in the case of pure iron injection at the source with a power law 
injection index y g = 2.2. Full squares correspond to the SimProp result while continuous lines correspond to the 
solution of the nuclei kinetic equation. Left panel: fluxes of iron and of secondary nuclei with A = 50,40. Right 
panel: fluxes of secondary nuclei with A = 30, 20, 10. 



10% with respect to the flux computed taking into account the intrinsic stochasticity of the process 

cm. 

In the case of nuclei propagation the energy losses due to the process of pair production can be 
simply related to the corresponding quantity for protons 1 4 1. Therefore the CEL approximation can be 
suitably used also for the pair production process involving nuclei. 

Let us now discuss the process of photo-disintegration of nuclei: this interaction changes the nu- 
cleus kind leaving its Lorentz factor unchanged. In the kinetic approach of H the process of photo- 
disintegration is treated as a decay process that simply depletes the flux of the nucleus A. The photo- 
disintegration "life-time" is defined as: 

with A the atomic mass number and r the Lorentz factor of the interacting particle, e' the energy of 
the background photon in the rest frame of the particle, €o(A) the threshold of the considered reaction 
in the rest frame of the nucleus A, <x the relative cross section, v the photo-disintegration multiplicity 
(number of emitted nucleons), e the energy of the photon in the laboratory system and n y (e) the density 
of the background photons per unit energy. 

The kinetic approach is based on the transport equations that describe the propagation of protons 
and nuclei through astrophysical backgrounds, these equations for protons and nuclei are respectively: 

dnJr, t) d r n 

dT^ ~ dr [^ (r ' t)np{r ' °1 = Qp(r ' (2) 

dn A (r, t) d n A (r, t) 

of or T A (r,t) 

where n is the equilibrium distribution of particles, b are the energy losses (adiabatic expansion of 
the Universe and pair/photo-pion production for protons or only pair-production for nuclei) Q is the 
injection of freshly accelerated particles and, in the case of nuclei, also the injection of secondary 
particles produced by photo-disintegration (secondary nuclei). 

Concerning sources, in this paper, we will consider the simplest hypothesis of a uniform distri- 
bution of sources that inject UHECR of different species with a power law injection of the type: 
Qinj(r,z) = (2o(z)*r r9 <?~ r / r " rat , where r is the Lorentz factor of the particle, r max is the maximum 
Lorentz factor that the sources can provide and y CJ is the power law index of the distribution of the 
accelerated particles at the sources, being both parameters r max and y g independent of the particle 
type (nucleus or proton). The injection of secondary nuclei should be modeled taking into account 
the characteristics of the photo-disintegration process. The dominant process of photo-disintegration 
is one nucleon (N) emission, namely the process (A + 1) + jbkg —* A + N, this follows directly from 
the behavior of the photo-disintegration cross-section (see [4| and references therein). Moreover, at 
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Fig. 2. UHECR flux compared with HiRes and Auger data. Left panel: dip model and HiRes data, pure proton 
injection with a power law index at the source y 9 = 2.7 and a maximum acceleration energy E p max = 10 21 eV. 
Right panel: disappointing model and Auger data, maximum proton energy = 4 x 10 18 eV with a mixture of 
p, helium and iron at the source, the injection power law index is y g = 2.0. 



the typical energies of UHECR (E > 10 17 eV) one can safely neglect the nucleus recoil so that 
photo-disintegration will conserve the Lorentz factor of the particles. Therefore the production rate 
of secondary A-nucleus and A-associated nucleons will be given by Qa(E,z) = Qp(r,z) = 
where Ta+i is the photo-disintegration life-time of the nucleus father (A + 1) and n^+i is its equilibrium 
distribution, solution of the kinetic equation for the nucleus A + 1 . 

The kinetic equations (12131 1 together with the injection of primary and secondary nuclei can be 
solved analytically obtaining the equilibrium distribution of primary injected nuclei and all kind of 
secondaries produced along the photo-disintegration chain. 

Let us now discuss the implementation of the stochastic treatment of the nuclei propagation as 
developed in the SimProp MC approach IfTTI . The calculations of the energy evolution and of the 
photo-disintegration life-time (Eq. ([1)) are performed step by step in red-shift. The survival probability 
as a function of red-shift and Lorentz factor of the nucleus A is: 

P(r,z) = expf- P —I— 

where z and z* are the values of the redshift of the current step (from z* to z). In the standard cosmology 
the term \dt/dz\ is given by|^| = (1+ ' zx , H(z) - Hq^(1 + z) 3 £2 m + &a with Ho,Q m ,QA chosen 
according to the standard cosmology Ifl2ll . 

In the SimProp MC code only the photo-disintegration process is treated stochastically while the 
other processes of energy losses are treated as in the kinetic approach, i.e. under the CEL approxima- 
tion. 

Figure[T]shows the good agreement of the fluxes computed with the MC approach and in the kinetic 
approach, for this comparison we have chosen an ideal case of a pure Fe injection at the source with 
a power law injection index y g - 2.2. In Ml 1141 a more detailed comparison of the results obtained 
with different computation schemes was performed, three different MC schemes were compared each 
other and with the kinetic approach. The fairly good agreement of the fluxes computed with different 
computation schemes shows the solidity of the theoretical approaches used to describe the propagation 
of UHE particles, giving a very useful workbench to test different theoretical models for UHECR 
sources. 



at 
dz 



dz 



(4) 



3 Transition from Galactic to Extra-Galactic Cosmic Rays 

In this section we will consider two different models aiming to explain the experimental observations 
of UHECR. In the case of HiRes and TA observations the experimental data are consistent with a pure 
proton composition, this experiments favor the dip model that is based on the hypothesis of a pure 
proton injection at the sources with an injection power law index y g = 2.7 [7|. The Auger observations 
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Fig. 3. Transition from galactic to extra-galactic cosmic rays. The galactic cosmic ray flux is taken from 1151 . 
The experimental data (as labeled) are those form HiRes-MIA |10|, Auger |10|, Kaskade [16| and the average 
over different experimental results as in JT7 ]. Fluxes of the extragalactic components are obtained with the same 
assumptions of figure|2] Left panel: dip model. Right panel: disappointing model. 



favor a mixed composition at the highest energies with a substantial nuclei pollution in the spectrum 
already at energies E > 4 x 10 18 eV. The Auger observations on fluxes can be accommodated in the so- 
called disappointing model 11131 . that is based on the assumption of a flat injection y g = 2.0 and a mixed 
composition at the sources made of protons, helium and iron with an high energy suppression of the 
fluxes due to the maximum attainable injection energy. In the disappointing model primary injected 
protons are assumed with a relatively low maximum energy E^ mx — 4 x 10 18 eV and, in a rigidity 
dependent scenario, primary injected nuclei will have a maximum energy given by E^ WK = ZE^ mx , 
being Z the atomic mass number of the nucleus A. Therefore, in the disappointing model the behavior 
of the fluxes at the highest energies is fixed by the maximum energy and it is not a propagation feature 
(as the GZK suppression) lfl3l . In figure [2] we plot the comparison of the dip model with the HiRes 
data 1 10 1 (left panel) and of the disappointing model with the Auger data 1 10] (right panel). In the latter 
comparison we have included in the Fe curve the sum of primary iron and all secondaries with A > 40, 
while the sum of the remaining secondaries with A < 40 are plotted in the magenta curve below Fe. 
Figure [2] shows the agreement of the dip model with the HiRes observations and of the disappointing 
model with the Auger observations. This fact signals once more the puzzling situation we are facing 
in UHECR physics with different experiments favoring alternative scenarios. 

In order to have a deeper view of the problem, trying to asses the intimate nature of UHECR, 
it is useful to consider the transition among galactic and extra-galactic cosmic rays, comparing the 
matching of the theoretical fluxes with experimental data fl4"l . In the analysis that follows we will 
assume a single fixed model for the galactic component, namely the flux computed in lfT31l taking into 
account the large scale distribution of supernova remnants in the galaxy. In figure [3] we plot the total 
flux with the galactic component of lfT31l and the extra-galactic components obtained in the two cases 
of the dip model (left panel) and of the disappointing model (right panel). The experimental data of 
figure[3]are those of HiRes-MIA IflOll . Auger BIOL Kaskade [16] and the average of the fluxes measured 
by different experiments as presented in fT71 . The matching of galactic and extragalactic fluxes in the 
case of the dip model (figure [3] left panel) gives a very good description of the experimental data, 
reproducing in an extremely accurate way the spectra observations in the intermediate energetic regime 
where the transition stands. The case of the disappointing model gives a less accurate description of 
the experimental data (figure [3] right panel) with a slightly suppression of the theoretical flux in the 
transition region not seen experimentally. 



4 Conclusions 

The theoretical study of the propagation of UHECR through astrophysical backgrounds has reached 
a remarkable level of refinement with several different approaches providing very reliable results. 
Here, in particular, we have compared the kinetic approach of |4] with the MC approach presented 
in IfTTII showing the very good agreement of their results (figure[TJ. Using these theoretical tools it is 
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possible to study with good accuracy different models for the production of UHECR by comparing 
experimental data with theoretical expectations. 

In the present paper we have considered the dip model [7| and the disappointing model fT3"l . The 
dip model is based on the assumption of a pure proton composition of UHECR with an injection 
power law index y g > 2.5, the experimental data of HiRes and TA confirm with very good accuracy 
this model (figure [2j left panel) while the Auger data disfavor it with an heavier composition at the 
highest energies. The disappointing model, on the other hand, is based on the assumption of a mixed 
composition at the sources with an injection power law index y g < 2.5 and a relatively low maximum 
energy for protons (E' t ' mx = 4 x 10 18 eV). This model provides a very accurate description of the 
Auger data (figure |2l right panel), while it does not reproduce the HiRes and TA observations. To 
better characterize these theoretical models, we have also focused the attention on the transition among 
galactic and extra-galactic CR assuming the galactic component as computed in ifTBll . In the case of the 
dip model the transition is a sharp change from heavy (galactic) to light (extra-galactic) dominance, 
starting at energies around 10 17 eV. The total flux in the case of the dip model shows a very good 
agreement with the all particle spectrum observed by different experiments (figure[3| left panel). In the 
case of the disappointing model the transition is placed at energies around 10 18 eV and the total flux 
seems not in perfect agreement with observations in particular at energies around the transition region 
(figure [3] right panel). 
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